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Background: A subset of patients with chronic myeloid leukemia (CML) can sustain a complete 
molecular response after discontinuing imatinib mesylate (IM). We focused on microRNAs 
(miRNAs), with the aim of finding a molecular biomarker to discriminate which patients can 
safely and successfully discontinue IM use. 

Methods: To identify miRNAs that showed altered expression in patients who had discontinued 
IM (STOP-IM group), we first screened miRNA expression of peripheral blood mononuclear 
cells by using a TaqMan miRNA array on samples from five unselected patients from the STOP- 
IM group, seven CML patients receiving IM (IM group), and five healthy volunteers. We then 
performed miRNA quantification in 49 CML patients with deep molecular response. Mann- 
Whitney U and chi-square tests were used to determine statistical significance for comparisons 
between the control (healthy volunteers) and test groups (STOP-IM and IM groups). Multiple 
groups were compared by one-way analysis of variance. 

Results: Downregulation of miR-148b was noted in patients in the STOP-IM group and 
in a subset of the IM group. We then subdivided the IM patients into two groups: one with 
downregulated miR-148b expression (IM-1; less than the cut-off value) and the other without 
downregulated miR-148b expression (IM-2; greater than the cut-off value). The number of 
patients who had a sustained stable molecular response was significantly lower in IM-2 group. 
This group also had a significantly lower percentage of natural killer cells. 
Conclusion: Downregulated miR-148 may contribute to immune surveillance in STOP-IM 
patients and may therefore have potential as additive information in managing CML patients 
undergoing treatment with IM. 
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Background 

The treatment of chronic myeloid leukemia (CML) has been greatly impacted by the 
development of imatinib mesylate (IM), the first specific inhibitor of BCR-ABL1. 1 
In newly diagnosed patients with chronic-phase CML, IM induces a complete hema- 
tological response in more than 90% of cases and a complete cytogenetic response 
in more than 80% of cases, 1 and IM has been established as the standard therapy 
for CML. Recent reports have demonstrated that a subset of CML patients with a 
complete molecular response (currently referred to as undetectable minimal residual 
disease, UMRD) can sustain this state after IM discontinuation (STOP-IM). 2 ^ 1 We and 
others have postulated that immunological surveillance may play an important role in 
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these patients. 5-8 More recently, Ng et al demonstrated that 
East Asian CML patients with a common intronic deletion 
polymorphism in the gene coding BIM (at intron 2) had an 
inferior response to tyrosine kinase inhibitors compared to 
those without the genetic variation. 910 Taken together, immu- 
nological diversity as well as genetic background should 
be considered in discussions of the efficacy and safety of 
IM cessation in CML patients. Currently there is no defini- 
tive classifier to discriminate CML patients with UMRD who 
could safely stop IM and those who should not. 

MicroRNAs (miRNAs), which are noncoding single- 
stranded RNAs of 21-25 nucleotides, have recently been 
implicated in the regulation of cellular processes such as 
development, differentiation, apoptosis, proliferation, and 
hematopoiesis. MiRNAs have also been implicated in the 
development of human cancers, either as tumor suppressors 
or as oncogenes. 11 12 Approximately 100 miRNAs are known 
to be expressed in cells of the immune system and to play key 
roles in the development of innate and acquired immunity, 1314 
therefore, aberrant expression of miRNAs that target the 
immune system may be related to the immunological back- 
ground of cancer patients. In CML cells, abnormal expression 
of several miRNAs (miR-15a, miR-16, miR-31, miR-142, 
miR-155, miR-181, miR-221, miR-564, let7a, and the 
polycistronic miR-17-92 cluster) has been described. 1518 

Although the miRNA expression profile in CML cells 
has been extensively studied little is known about miRNA 
expression in CML patients with UMRD whose CML cells 
cannot be detected. We therefore determined the miRNA 
profiles of peripheral blood mononuclear cells (PBMCs) in 
CML patients with UMRD who were receiving IM or had 
discontinued IM, with the goal of clarifying whether altered 
miRNA expression is linked to pathophysiology in CML 
patients who safely discontinued IM. Here, we provide a 
distinct miRNA expression profile of circulating PBMCs in 
CML patients with successful IM cessation. We additionally 
describe a subset of CML patients receiving IM that showed 
similar miRNA expression profiles. 

Materials and methods 

Patients 

Forty-nine consecutive patients with CML who had UMRD 
were enrolled in this study. In 1 6 patients who had discon- 
tinued IM and had sustained UMRD for more than 6 months 
(STOP-IM group), samples were collected when IM was 
stopped. In contrast, 33 patients who were receiving IM and 
had UMRD at the time of sample collection were categorized 
as the IM group. Fifteen healthy volunteers served 



as controls. CML patients before treatment and those who 
did not respond to IM were not included in this study. There 
was no significant difference in sex or age between the CML 
groups and the healthy volunteers. UMRD was defined as the 
absence of a BCR-ABL transcript as determined by nested 
reverse transcription-polymerase chain reaction (RT-PCR) 
or a highly sensitive transcription-mediated amplification 
(TMA) method. These PCR methods can detect at least a 
4-log reduction in the BCR-ABL transcript (international 
scale <0.01%). The study was approved by the institutional 
review board of Tokyo Medical University (no 930, approved 
24 June 2008). Written informed consent was obtained from 
all participants prior to specimen collection, in accordance 
with the Declaration of Helsinki. 

RNA isolation from peripheral blood 

PBMCs were obtained from 1 0 mL of whole blood by the Ficoll 
density gradient method (GE Healthcare Bio-Sciences AB, 
Uppsala, Sweden). Cell pellets were quickly frozen at -80°C 
until use. Total RNA was isolated from cell pellets with a 
mirVana PARIS kit (Ambion, Austin, TX, USA), as reported 
previously. 19 Briefly, 600 uL of cell lysate was mixed with 
an equal volume of denaturing solution, and the subsequent 
phenol extraction and cartridge filtration were carried out 
according to the manufacturer's instructions. 

TaqMan low-density array screening 

The miRNA expression profiles were determined with a 
TaqMan Array Human microRNA Card A (Thermo Fisher 
Scientific, Waltham, MA, USA). Ninety nanograms of total 
RNA was used as an input in each RT reaction. The RT 
reaction and pre-amplification step was set up according to 
the manufacturer's recommendations. miRNAs were reverse 
transcribed with the Megaplex Primer Pools (Human Pools 
A version 2.1; Thermo Fisher Scientific). RT reaction prod- 
ucts from the PBMC samples were further amplified with 
Megaplex Pre Amp primers (Primers A version 2.1; Thermo 
Fisher Scientific). Quantitative RT-PCR was performed on 
an Applied Biosystems 7900HT thermocycler according to 
the manufacturer's recommended program. With the use of 
SDS 2.4 software and DataAssist (Thermo Fisher Scientific), 
the expression of miRNA was calculated based on cycle 
threshold (Ct) values normalized by those of RNU6B. Data 
analysis was done using GeneSifter* software (Perkin Elmer, 
Waltham, MA, USA). P-values of less than 0.05 were consid- 
ered to indicate a statistically significant difference, and the 
Benjamini-Hochberg algorithm was used for estimation of 
false discovery rates, as we have reported previously. 19 
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Real-time quantitative RT-PCR 

We determined the amounts of the individual miRNAs by real- 
time quantitative RT-PCR with a TaqMan MicroRNA Assay 
(Thermo Fisher Scientific) and the following miRNA-specific 
stem-loop primers: hsa-miR-148b (000471), hsa-miR-326 
(000542), hsa-let7b (000378), and RNU6B (001093). 
Subsequently, quantitative real-time PCR was performed 
with an ABI Prism 7000 sequence detection system (Thermo 
Fisher Scientific). The reaction was initiated by incubation 
at 95°C for 2 minutes, followed by 50 cycles of 95°C for 
1 5 seconds and then 60°C for 1 minute. All reactions were run 
in duplicate. Mean (Ct) values for all miRNAs were quanti- 
fied with sequence detection system software (SDS version 
1.02; Thermo Fisher Scientific). All miRNA expression was 
normalized to RNU6B expression, yielding a -ACt value, as 
reported elsewhere. 20 The -AACt value was then calculated 
by subtracting the -ACt value of the normal sample from the 
respective -ACt values of patient samples. Expression of all 
miRNAs was normalized by using the 2~ AACt method. 

Flow cytometry 

Freshly obtained PBMCs were separated using a Ficoll den- 
sity gradient (GE Healthcare Bio-Sciences) and resuspended 
in Roswell Park Memorial Institute (RPMI) 1640 medium 
supplemented with 10% fetal calf serum. Immunophenotyping 
was performed on a FACSCanto II (BD Biosciences, San Jose, 
CA, USA) flow cytometer using CellQuest software (BD 
Biosciences) with a 5-color flow cytometry panel including 
antibodies against the following cell surface antigens and 
effector molecules: CCR7, CD3, CD8, CD45RO, and CD56 
(eBioscience, San Diego, CA, USA, and BioLegend, San 
Diego, CA, USA). Cells were analyzed with the FACSCanto II 
(BD Biosciences) flow cytometer using CellQuest software. 8 

Statistical analysis and pathway analysis 

Data are expressed as means + standard deviation (SD). 
Mann- Whitney U and chi-square tests were used to deter- 
mine statistical significance for comparisons between the 
control and test groups. Multiple groups were compared 
by one-way analysis of variance. GraphPad Prism software 
(version 5c for Macintosh; GraphPad Software Inc., La Jolla, 
CA, USA) was used for statistical analyses. Results were 
considered statistically significant when P<0.05. 

Following identification of differentially expressing 
miRNAs, the predicted targets for these altered miRNAs were 
subjected to micro RNA Target Filter in Ingenuity Pathways 
Analysis (IPA) software (Ingenuity System, Redwood City, 
CA, USA). 



Results 

The miRNA expression profiling 
byTaqman miRNA array 

To identify candidate miRNAs that showed altered expres- 
sion in patients in the STOP-IM group, we first screened 
miRNA expression using a TaqMan miRNA array (Thermo 
Fisher Scientific) on five unselected patients from the 
STOP-IM group, seven from the IM group, and five healthy 
volunteers. Among these three groups, we observed dif- 
ferential expression of 22 miRNAs as identified by using 
GeneSifter software (Figure 1 ; Gene Expression Omnibus 
[GEO] accession no GSE47652). Although 22 miRNAs were 
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Figure I miRNA profiling by TaqMan (Thermo Fisher Scientific, Waltham, MA, 
USA) miRNA array. 

Notes: A distinct expression pattern was found in the STOP-IM group (middle lane). 
Using SDS (version 2.4) and DataAssist software (Thermo Fisher), we calculated the 
expression of miRNAs based on their Ct values normalized by the Ct value of RNU6B. 
Data were analyzed with GeneSifter software {Perkin Elmer, Waltham, MA, USA). 
P-values of <0.05 were considered to indicate a statistically significant difference, and 
the Benjamini— Hochberg algorithm was used for estimation of false discovery rates. 
Abbreviations: miRNA, microRNA; Ct, cycle threshold; IM, imatinib mesylate; 
STOP-IM, patients who had discontinued IM and had sustained UMRD for more 
than 6 months; UMRD, undetectable minimal residual disease. 
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selected based on fold change, the expression profile was 
highly variable among individuals. We therefore selected 
three of these miRNAs (let-7b, miR-148b, and miR-326) 
for further validation on the basis of their averaged fold 
change in expression (log2 variation >0.6), significant 
P-value, quality control value, and their predicted target 
genes (Table SI). 

Quantification of individual miRNA 
by real-time quantitative RT-PCR 

We then validated the miRNA array results by real-time 
quantitative RT-PCR using samples from 16 patients of 
the STOP-IM group, 33 of the IM group, and 15 healthy 
volunteers. Validation of miR-148b, miR-362, and let-7b 
revealed that expression of miR-148b was significantly 
lower in the STOP-IM group than in the healthy volunteers 
(Figure 2A). The miR-148b expression tended to be lower 
in the STOP-IM group than in the IM group, although the 
difference was not statistically significant. We observed no 
statistically significant difference in the expression of let-7b 
between the STOP-IM and IM groups (Figure 2B). Despite 



miR-362 being upregulated in the STOP-IM group based 
on miRNA array screening, the validation analysis using a 
large number of samples revealed that miR-326 tended to be 
downregulated in the STOP-IM group in comparison to the 
IM group and the control group. We observed no statistically 
significant difference in the expression of miR-326 between 
each group (Figure 2C). 

Clinical and biological relevance 
of miR- 1 48b expression 

Based on the statistical significant results from quantifying 
individual miRNAs, we particularly focused on miR- 148b. 
To estimate a cut-off value for miR- 148b expression in 
mononuclear cells for distinguishing between the STOP-IM 
and healthy control groups, we performed receiver operat- 
ing characteristic (ROC) curve analysis (Figure 2D). The 
area under the ROC curve was 0.8792 with 75% sensitiv- 
ity (95% confidence interval [CI], 47.62%-90%) and 80% 
specificity (95% CI, 51.91%-95%). The cut-off value was 
thus determined to be 0.5880, which was used in subsequent 
analyses of miR- 1 48b expression. We next subdivided the IM 
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Figure 2 Expression of miR- 1 48b (A), let-7b (B), and miR-326b (C) in the STOP-IM, IM, and control groups. (D) Cut-off level for miR- 1 48b expression between the STOP- 
IM and control groups. 

Notes: The asterisk indicates a statistically significant difference. The area under the ROC curve was 0.8792 with 75% sensitivity (95% CI, 47.62%-90%) and 80% specificity 
(95% CI, 5 1 .9 1 %-95%). The cut-off value was thus determined to be 0.5880, which was used in the following analyses of miR- 1 48b expression. 

Abbreviations: IM, imatinib mesylate; ROC, receiver operating characteristic; CI, confidence interval; CML, chronic myeloid leukemia; STOP-IM, patients who had 
discontinued IM and had sustained UMRD for more than 6 months; UMRD, undetectable minimal residual disease. 
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patients into two groups: one with downregulated miR-148b 
expression (IM-1; less than the cut-off value) and the other 
without downregulated miR-148b expression (IM-2; greater 
than the cut-off value) (Table 1). 

Comparison of the STOP-IM, IM-1, and IM-2 groups 
revealed no significant differences in sex, age, Sokal score, 
duration of IM therapy, months to achieve UMRD from 
diagnosis, or prior interferon (IFN)-a therapy (Table 1 ). The 
total IM dose was significantly higher (P=0.0289) in the IM-2 
group (893+294 g) than in the STOP-IM group (636±369 g), 
while there was no statistically significant difference between 
the IM-1 group and STOP-IM group. Similarly, the number 
of patients who sustained stable molecular response was 
significantly lower (,P=0.0232) in the IM-2 group than in the 
STOP-IM group. Notably, the percentage of natural killer 
(NK) cells was significantly lower (P=0.03) in the IM-2 group 
than in the STOP-IM group. We could not find any statistical 
significance of BIM deletion among the 3 groups. 

Network analysis of miRNA- 148b 

ROC analysis of the miR-148b suggests that the miRNA 
signature in PBMCs may predict whether IM discontinuation 



Table I Clinico-hematological, immunological, and molecular 
characteristics of patients with and without downregulation of 
miR- 1 48b expression 



Downregulation 


STOP-IM 


IM-1 


IM-2 


P-value 


of miR-l48b 


Yes 


Yes 


No 




expression 










Number of patients 


16 


18 


15 




Age (years) 


57.5+1 1 .6 


60.4±9.l 


54.5±I5.4 


0.404 


Sex (male/female) 


7/9 


12/6 


1 1/4 


0.201 s 


Sokal score 


0.74±0.22 


0.70±0.09 


0.73±0.24 


0.849 


Total imatinib dose 


637±360 


848±324 


9541245 


0.0579 


(grams) 


1 


0.018511 
















Duration of imatinib 


67.4±34.9 


75.8±29.3 


88.9176.0 


0.2264 


therapy (months) 










Months to achieve 


S9.3±68.0 


48.2±39.8 


67.5144.9 


0.5473 


UMRD from diagnosis 










Stable MR* (yes/no) 


14/2 


13/5 


7/8 


0.0454* 




1 


0.023211 . 


1 




Prior IFN-a (yes/no) 


4/12 


3/15 


5/10 


0.3888* 


NK cells (%) 


I3.2±6.9 


I2.2±9.l 


8.110.96 


0.141 




1 


0.03H 


1 




CCR7-CR4SRO- 


8.4±3.9 


I0.5±5.4 


7.616.1 


0.1733 


CD8 + cells (%) 










BIM deletion (yes/no) 


1/ 10 


2/16 


1/13 


0.9288* 



Notes: *lf BCR-ABLl remained undetectable for a continuous period of 24 months 
of imatinib treatment, the response was termed stable molecular response (MR). 
Patients who received IM and did not have a sustained UMRD but maintained a 
major molecular response for more than 24 months (fluctuating UMRD) were 
not referred as "stable MR". P values are calculated by one-way ANOVA unless 
otherwise indicated; ^Student's t-test; 5 chi-square test. 

Abbreviations: IM, imatinib mesylate; UMRD, undetectable minimal residual disease; 
ANOVA, analysis of variance; IFN, interferon; NK, natural killer; STOP-IM, patients 
who had discontinued IM and had sustained UMRD for more than 6 months. 



is appropriate. Therefore, using in silico network analysis we 
assessed whether predicted targets of miR- 148b signature could 
provide insights into analogous changes enriched in the PBMC 
signaling pathways in CML-UMRD patients. Network analysis 
on predicted targets of miRNA- 148 was carried out using IPA 
(Ingenuity System), which provides an independent method 
of assessing the role of predicted miR-148b targets. The top 
network has ten nodal molecules that are predicted targets of 
miR-148b like DNMT3B, NR1I2, andTLR3 (Figure 3). 

Discussion 

One of the major questions in research on STOP-IM centers 
is the genetic and epigenetic markers for discriminating 
patients who can safely stop IM. Several empirical clini- 
cal studies have shown that durable stable UMRD before 
discontinuation, prior IFN therapy, and early deep molecular 
response are important to avoid molecular relapse. 2 - 7,21,22 In 
this context, we particularly focused on miRNA expression 
profiling of PBMCs in CML patients who could safely stop 
IM. Since all patients in this study exhibited deep molecular 
response, we did not do miRNA profiling of CML cells; there- 
fore, we could compare CML patients with or without IM and 
healthy volunteers. We found a distinct miRNA profile that 
was specific for the STOP-IM group, and we concluded that 
the miR- 148b expression level might be a possible diagnostic 
marker for safely stopping IM. 

miR-148a, miR-148b, and miR-152 are the three mem- 
bers of the miR148/152 family. 23 Aberrant expression of the 
miR148/152 family has been reported in tumor as well as 
non-tumor disease. The present study was performed using 
non-tumor PBMCs; target genes regulated by miR-148b 
might be different from those in tumor tissue. For this reason, 
we tried to identify possible targets through use of the IPA 
software (Ingenuity System). There were several candidate 
targets for miR- 14 8b as shown in Figure 3, such as chole- 
cystokinin B receptor (CCBR), DNA methyltransferase 1 
(DNMT1), DNMT3, nuclear receptor subfamily 1, group 1, 
member 2 (NR1 12), peroxisome proliferator-activated recep- 
tor alpha (PPARA), ribosomal protein S6 kinase, and poly- 
peptide 5 (RPS6KA). Notably, DNMT1 and DNMT3 were 
among the candidates. The CpG island in the first intron of 
Foxp3 DNA is demethylated in naturally occurring human 
and murine T-regulatory (Treg) cells, and the methylation 
status is linked with Foxp3 expression. 24 A recent report by 
Wang et al demonstrated Foxp3 + Treg cells require Dnmtl 
expression to prevent development of lethal autoimmunity. 23 
Therefore, we could not rule out the possibility that down- 
regulation of miR- 148b induces DNMT1 expression and 
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Figure 3 Network analysis of miR-l48b by Ingenuity Pathways Analysis. 

Notes: The miR-l48b-3p is located in the center. Solid lines indicate a direct relation, and broken tines indicates an indirect relation. 



promotes Treg suppression regulated by Foxp3 expres- 
sion. In the present study, we were not able to analyze Treg 
function; however, Rea et al recently reported that NK cells 
rather than Treg cells are associated with clinical outcome 
after IM discontinuation. 26 Since imatinib itself is known to 
inhibit CD4 + CD25 + Treg activity, 27 regulation of Treg after 
IM discontinuation might be partly compensated by down- 
regulation of miR-148b in a subset of STOP-IM patients. 

In addition to DNMT 1 , we noted that miR- 1 48b is indi- 
rectly linked to toll-like receptor 3 (TLR3) in the network 
presented by IPA (Ingenuity System). Liu et al recently 
reported that miRNA-148 and miRNA-152 fine-tune the 
regulation of the innate immune response and the antigen- 
presenting capacity of dendritic cells, and this fine-tuning 
may contribute to immune homeostasis and immune regu- 
lation. 28 The expression of miR- 148a, 148b, and miR152 
was upregulated in mouse dendritic cells upon maturation, 
and activation was induced by Tlr3, Tlr4, and Tlr9 agonists. 
In the clinical context, we subdivided CML patients based 
on miR- 148b expression and found that CML patients in 
the STOP-IM group and the IM-2 group showed different 
immunological features. Additionally, the percentage of NK 
cells was significantly lower in the IM-2 group. Although 
we did not intend to focus on immunological surveillance 
in STOP-IM patients at the beginning of this miRNA study, 
the present data may support previous reports regarding a 
possible role of NK cells in STOP-IM patients. 810 - 26 

In line with this, CML patients in the IM- 1 group could stop 
IM, and miR- 148 expression in PBMCs might be an additive 
biomarker for IM cessation. The correlation between miR- 1 48b 



levels and treatment- free remission is not clear at this moment 
(Figure S 1 A); therefore, follow-up study is needed to clarify the 
clinical significance of miR- 1 48b. We only had two patients who 
relapsed shortly after IM discontinuation (Figure SIB), and all 
STOP-IM patients enrolled in the current study had sustained 
UMRD for more than 12 months. Consequently, we could not 
determine whether miR- 1 48b levels are correlated with relapse 
after stopping IM. Since the number of patients in the present 
study was rather small, an analysis of a larger number of patients 
in an independent cohort may clarify this point. 

Furthermore, it is a matter for discussion whether PBMCs 
are suitable for miRNA analysis in clinical specimens. Analysis 
of the lymphocyte fractions (eg, NK cells and T cells) would be 
useful, but using fractionated samples from clinical specimens 
is difficult in practice. For this reason, we used circulating 
lymphocytes rather than purified NK or T cells; however, we 
could not completely rule out the possibility that the differ- 
ence in miRNA expression was due to differences in cellular 
composition. Another limitation in the current study is that 
the miRNA profile was not from cancer cells; therefore, we 
could not perform a functional analysis of the target genes. 
Unfortunately, we could not perform biological experiments 
to demonstrate a possible association between downregula- 
tion of miR- 148b and the immune surveillance in CML. We 
alternatively applied an IPA (Ingenuity System) bioinformat- 
ics tool to estimate the network between miR- 148b and its 
target genes (Figure 3). Although our bioinformatics analysis 
is mainly based on previous knowledge base, it may help to 
understand the biological relevance of miR- 148b in the cur- 
rent study. Finally, we could not find any relationship between 
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BIM deletion and the level of miR- 148b expression. This point 
should be clarified in a larger cohort in the near future. 

Conclusion 

Our observation of downregulated expression of miR- 148b 
supports the idea that miRNA may contribute to immune sur- 
veillance in STOP-IM patients and may therefore have potential 
as an additive information marker in managing CML patients 
being treated with IM. Whether miR-148b may be used as a 
biomarker among other factors (eg, NK cells) to predict disease 
relapse after IM discontinuation needs to be investigated. 
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TFS (months) Months 

Figure SI (A) miR-l48b expression and treatment-free survival (TFS) in STOP-IM patients. Two patients who relapsed shortly after IM discontinuation are shown in open 
circles. (B) Sequential analysis of miR-l48b expression in two patients who relapsed shortly after IM discontinuation. 

Note: Patients who had discontinued IM and had sustained UMRD for more than 6 months (STOP-IM group). In contrast, patients who were receiving IM and had UMRD 

at the time of sample collection were categorized as the IM group. 

Abbreviation: IM, imatinib mesylate; UMRD, undetectable minimal residual disease. 
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0.027347 
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Table S I Fold changes, P-values, and quality control values of 22 genes differentially expressed among 3 groups (GEO accession number GSE47652). 
Abbreviations: ANOVA, analysis of variance ID, identification document; SEM, standard error of measurement; GEO, Gene Expression Omnibus. 
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